Axial charges of N(1535) and N(1650) in lattice QCD with two flavors of dynamical 

quarks 
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We show the first lattice QCD results on the axial charge g 1 ^ of iV*(1535) and A r *(1650). The 
measurements are performed with two flavors of dynamical quarks employing the renormalization- 
group improved gauge action at /3=1.95 and the mean- field improved clover quark action with 
the hopping parameters, k=0.1375, 0.1390 and 0.f400. In order to properly separate signals of 
A r *(1535) and JV*(1650), we construct 2x2 correlation matrices and diagonalize them. Wraparound 
contributions in the correlator, which can be another source of signal contaminations, are eliminated 
by imposing the Dirichlet boundary condition in the temporal direction. We find that the axial 
charge of iV*(1535) takes small values as ~ O(0. 1), whereas that of iV*(1650) is about 

0.5, which is found independent of quark masses and consistent with the predictions by the naive 
nonrelativistic quark model. 



Introduction. Chiral symmetry is an approximate 
global symmetry in QCD, the fundamental theory of the 
strong interaction; this symmetry together with its spon- 
taneous breaking has been one of the key ingredients in 
the low-energy hadron or nuclear physics. Due to its 
spontaneous breaking, up and down quarks, whose cur- 
rent masses are of the order of a few MeV, acquire the 
large constituent masses of a few hundreds MeV, and 
are consequently responsible for about 99% of mass of 
the nucleon and hence that of our world. Thus one 
would say that chiral condensate {tpip}, the order pa- 
rameter of the chiral phase transition, plays an essential 
role in the hadron-mass genesis in the light quark sec- 
tor. On the other hand, chiral symmetry gets restored in 
systems where hard external energy scales such as high- 
momentum transfer, temperature(T), baryon density and 
so on exist, owing to the asymptotic freedom of QCD. 
Then, are all hadronic modes massless in such systems? 
Can hadrons be massive even without non- vanishing chi- 
ral condensate? 

An interesting possibility was suggested some years ago 
by DeTar and Kunihiro who showed that nucleons 
can be massive even without the help of chiral conden- 
sate due to the possible chirally invariant mass terms, 
which give degenerated finite masses to the members in 
the chiral multiplet (a nucleon and its parity partner) 
even when chiral condensate is set to zero: To show this 
for a finite-T case, they introduced a linear sigma model 
which offers a nontrivial chiral structure in the baryon 
sector and a mass-generation mechanism completely and 
essentially different from that by the spontaneous chiral 
symmetry breaking, fntcrestingly enough, their chiral 
doublet model has recently become a source of debate as 
a possible scenario of observed parity doubling in excited 
baryons 0, S H H , although their original work [l| 
was supposed to be applied to finite-T systems. 

It is thus an intriguing problem to reveal the chiral 
structure of excited baryons in the light quark sector be- 



yond model considerations. One of the key observables 
which are sensitive to the chiral structure of the baryon 
sector is axial charges The axial charge of a nucleon 
N is encoded in the three-point function 

(N\A°\N) = UY^59A(q 2 ) + q^l 5 h A (q 2 )]u. (1) 

Here, = Qj^js^Q is the isovector axial current. The 
axial charge <?a is defined by <?a(9 2 ) with the vanishing 
transferred momentum q 2 = 0. It is a celebrated fact that 
the axial charge g% N of 7V(940) is 1.26. Though the axial 
charges in the chiral broken phase can be freely adjusted 
with higher-dimensional possible terms and cannot be 
the crucial clues for the chiral structure [H, 0] , they would 
surely reflect the internal structure of baryons and would 
play an important role in the clarification of the low- 
energy hadron dynamics. 

In this paper, we show the first unquenched lattice 
QCD study [8] of the axial charge g%' N * of TV* (1535) and 
TV* (1650). We employ 16 3 x 32 lattice with two flavors 
of dynamical quarks, generated by CP-PACS collabora- 
tion Q with the renormalization-group improved gauge 
action and the mean-field improved clover quark action. 
We choose the gauge configurations at j3 = 1.95 with the 
clover coefficient csw = 1-530, whose lattice spacing a is 
determined as 0.1555(17) fm. We perform measurements 
with 590, 680, and 680 gauge configurations with three 
different hopping parameters for sea and valence quarks, 
K scaj K vai = 0.1375,0.1390 and 0.1400, which correspond 
to quark masses of ~ 150, 100, 65 MeV and the related 
ir-p mass ratios are mps/my = 0.804(1), 0.752(1) and 
0.690(1), respectively. Statistical errors are estimated by 
the jackknife method with the bin size of 10 configura- 
tions. 

Our main concern is the axial charges of the negative- 
parity nucleon resonances iV*(1535) and iV*(1650) in 
^ channel. We then have to construct an optimal 
operator which dominantly couples to iV*(1535) or 
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TV* (1650). We employ the following two independent 
nucleon fields, Ni(x) = e a h c u a {x){u b {x)C"f^d c {x)) and 
N 2 (x) = e a t lc '-f 5 u a (x)(u b (x)Cd c (x)), in order to construct 
correlation matrices and to separate signals of iV*(1535) 
and AT*(1650). (Here, u(x) and d(x) are Dirac spinor for 
u- and d- quark, respectively, and a, 6, c denote the color 
indices.) Even after the successful signal separations, 
there still remain several signal contaminations mainly 
because lattices employed in actual calculations are finite 
systems: Signal contaminations (a) by scattering states, 
(b) by wraparound effects. 

Comment to (a) : Since our gauge configurations 
are unquenched ones, the negative parity nucleon states 
could decay to it and N, and their scattering states could 
come into the spectrum. The sum of the pion mass M n 
and the nucleon mass Mjv is however in our setups heav- 
ier than the masses of the lowest two states (would-be 
A 7 "* (1535) and AT* (1650)) in the negative parity channel. 
We then do not suffer from any scattering-state signals. 

Comment to (b) : The other possible contamina- 
tion is wraparound effects [10(. Let us consider a two- 
point baryonic correlator (AT* (t sn k)N* (£ src )) in a Eu- 
clidean space-time. Here, the operators N*(t) and 
N*(t) have nonzero matrix elements, (0\N*(t)\N*) and 
(N*\N*(t)\0), and couple to the state \N*). Since we 
perform unquenched calculations, the excited nucleon N* 
can decay into N and 7r, and even when we have no scat- 
tering state \N + 7r), we could have another "scattering 
states". The correlator (N*(t sn k)N*(t SIC )) can still ac- 
commodate, for example, the following term. 



(Tr\N*(t snk )\N)(N\N*(t STC )\7r) 

x e -Bjv(* s „k-t S ro) x e -~E„(N t -t aak +t B 



(2) 



Here, N t denotes the temporal extent of a lattice. Such 
a term is quite problematic and mimic a fake plateau 
at En — E n in the effective mass plot because it be- 
haves as ~ e -( B N--E T )(« S nk-t S rc)^ Although these con- 
taminations disappear when one employ enough large- 
st lattice, our lattices do not have so large Nt. In 
order to eliminate such contributions, we impose the 
Dirichlct condition on the temporal boundary for valence 
quarks, which prevents valence quarks from going over 
the boundary. Though the boundary is still transpar- 
ent for the states with the same quantum numbers as 
vacuum, e.g. glueballs, such contributions will be sup- 
pressed by the factor of e ~ E oN t an( j we ne gl e ct them in 
this paper. (Wraparound effects can be found even in 
quenched calculations (Io||.) 

Formulation. We here give a brief introduction to 
our formulation HHIIi- Let us assume that we have a 
set of N independent operators, 0^ nk for sinks and 0^ c 
for sources. We can then construct an N x N correla- 
tion matrix C IJ (T) = (0^)0^(0)) = C , s t nk A(T)C src . 
Here, (Cj nk )« = (0|Of nk |») and (C^j = {j\0& |0> are 
general matrices, and A(T)y is a diagonal matrix given 



by K{T)ij = 5ije~ EiT . The optimal source and sink op- 
erators, Oll c and C* nk , which couple dominantly (solely 
in the ideal case) to i-th lowest state, are obtained as 
OH = EjO&(C m )J and OU = E;(CL)tf °s*U, 
since (C s t nk )- 1 C(T)(C* src )- 1 = A(T) is diagonal. Besides 
overall constants, (C src ) _1 and (C^) -1 are obtained as 
the right and left eigenvectors of C _1 (T + l)C(T) and 
C(T)C{T+ I)" 1 , respectively. 

The zero-momentum-projected point-type operators, 



Ni(t) 



£abcW Q (x, i)( w b (x, i)C 75 cf(x, t)) (3) 



and 



N 2 (t) 



£abc75« a (x, i)(u b (x, i)C<f (x, t)), (4) 



are chosen for the sinks. For the sources, we employ the 
following wall-type operators in the Coulomb gauge, 



AM*) = 



and 



xi,x 2 ,x 3 



£ ab cM Q (xi, i)(u 6 (x 2 , t)C lb J C (X 3 , t)) 



(5) 



N 2 (t)= ^ e abc75 S a (x 1 ,t)(M 6 (x 2 ,t)Cd c (x 3 ,t)). 

xi,x 2 ,x 3 

(6) 

The parity is flipped by multiplying the operator by 75; 
N+(t) = Ni(t) and N~(t) = 7 5 AT i (*). The optimized sink 
(source) operators A/^ (A/^ ± ), which couple dominantly 
to the i-th lowest state are constructed as 

A/?(t) = Nt(t)+[(C^/(C^]Nt(t) (7) 

= N?{t)+LfN±(t), (8) 

and 

N±{t) =Nt(t)+[(Ctro)2l/(Ctro)u]W(t) (9) 

= Nf(t) + RfNt(t). (10) 

Now that we have constructed optimized operators, 
we can easily compute the (non-renormalized) vector 
and axial charges gy^f*' for the positive- and negative- 
parity nucleons via three-point functions with the so- 
called sequential-source method [13] ■ In practice, we 
evaluate gy A (*) defined as 



±[lat] (;) = 



9v.A 



Tr Ty A (B(t snk )J^ A {t)B{t SIC )) 
Tr r v , A (B(t SDk )B(t sl . c )) 



(11) 



and extract gy^ by the fit gy^ = in the 

plateau region. B{t) denotes the (optimized) interpolat- 
ing field for nucleons, and Ty.A are 7 AI : ^r i and 7^5^^-, 
respectively. JY' A {t) are the vector and the axial vector 



3 



S -4 




o ok=0.1375 

B a K=0.1390 

o o k=0.1400 



2 4 6 

current insertion point t 



FIG. 1: The non-renormalized axial charge of iV*(1535), 
<;^ ' lat ' (t), as a function of the current insertion time t. 

currents inserted at t. We show in Fig. Q] g A °^ at ^ (t) for 
iV*(1535) as a function of the current insertion time t. 
They are rather stable around t src < t < i sn k- 

We finally reach the renormalized charges 
with the prefactors Z A y = 

2kuqZv,a (1 + i>v,A^ J ! which are estimated with the 
values listed in Ref. 0. 

We show the fitted values of Lf 2 and Rf 2 in Table |U 
L(T) and R(T) are rather stable and show a plateau 
from relatively small value of T (T ~ 2), which is the 
same tendency as that found in Ref. We plot in 



^ ±[lat] 

9a,v = z A,vg A ,v 



Fig. [2] Lf (T) and Rf(T) obtained at k=0.1390, for the 
purpose of reference. The energies E^ 2 are extracted 
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FIG. 2: As typical examples, L 1 and R 1 obtained at 
k=0.1390 are plotted. 

from two-point correlation functions by the exponential 
fit as (Mt(t sxc + T)jjf(t SIC )) = CeM-EfT) in the 
large- T region. The value at each hopping parameter 
is found to coincide with that in the original paper by 
the CP-PACS collaboration 0], with deviations of 0.1% 
to 1%. We here perform simple linear chiral extrapola- 
tions for E^ 2 . The chirally extrapolated values as well 
as those at each hopping parameter for E^ 2 in lattice 



unit are listed in Table |TJ Although the mass E+ of the 
ground-state positive-parity nucleon at the chiral limit is 
overestimated in our analysis (a -1 = 1.267 GeV), this 
failure comes from our simple linear fit. 
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FIG. 3: The renormalized vector and axial charges of the 
positive- and the negative-parity nucleons are plotted as a 
function of the squared pion mass m\. upper panel: The 
results of the vector charges. The solid line is drawn at gv — 1 
for reference, lower panel: The results of the axial charges. 
The solid line is drawn at gA = 1-26 and the dashed line is 
drawn at gA = 0. 

Results. We first take a stock of the vector charges 
of the ground-state positive- and negative-parity nu- 
cleons as well as the axial charge g A + of the ground-state 
positive-parity nucleon, which are well known and can 
be the references. We show gy , the vector charges of 
the positive- and the negative-parity nucleons obtained 
with three hopping parameters k=0.1375, 0.1390 and 
0.1400, in the upper panel in Fig. [31 where the verti- 
cal axis denotes g^ 1 and the horizontal one the squared 
pion masses. (These values are also listed in Table HI) 
The vector charges should be unity if the charge conser- 
vation is exact, whereas we can actually find about 10% 
deviations in Table U or in the upper panel in Fig. [3] 
Such unwanted deviations are considered to arise due to 
the discretization errors: The present lattice spacing is 
about 0.15 fm, which is far from the continuum limit. In 
fact, the decay constants obtained with the same setup 
as ours deviate from the continuum values by 0(10)%. 
We should then count at least 10% ambiguities in our re- 
sults. The axial charge g A + of the positive parity nucleon 
is also shown in the lower panel in Fig.[3J As found in the 
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TABLE I: upper table: The fitted values of L* 2 and Rf 2 f° r the ground and the 1st excited states in positive- and 
negative-parity channels are listed. M n and E^ 2 denote the pion mass and the ground- and the 1st excited-state energies for 
the positive- and the negative-parity channels for each k in the lattice unit, (a = 0.1555fm and a -1 = 1.267GeV.) The row 
"C.L." shows the values at the chiral limit, lower table: The non-renormalized vector and axial charges for n-th positive- 
and negative-parity nucleons ^™^( ,1 ' d )[ lat ] are ]i s ted. The superscripts (u) and (d) denote the u- and d-quark contributions, 



respectively. The total axial and vector charges, gy^' 13 *' 
as the renormalization factor and the improvement coefficients Zv.a 



n±(u)[lat] 
9v 



n±(d)[lat] 

9v 



and g A 
2kuoZv,a (l + bv,A 



±[lat] 



are also listed, J 



Zv,A9va denote the renormalized charges. 



n±(n)[lat] 
9 A 



s n±(d)[l a t] ) agwell 



for vector and axial currents 



Rt 



R; 



EI 



El 



e; 



0.1375 -0.4341 -0.4573 0.0355 0.0126 -1353 -314.1 -1.432 
0.1390 -0.4526 -0.4552 0.1115 -0.2036 -845.9 -228.1 -2.729 
0.1400 -0.1605 -0.3552 0.0990 -0.0151 -408.9 -143.6 -1.510 



-1.302 0.8985(5) 1.696(1) 2.137(10) 2.524(53) 2.141(14) 
-1.084 0.7351(5) 1.459(1) 1.854(13) 2.162(44) 1.908(17) 
-1.038 0.6024(6) 1.270(2) 1.665(15) 2.046(67) 1.733(25) 



C.L. 



0.936(3) 1.277(25) 1.570(109) 1.411(38) 



0+Tlat] 

9y 



0— [latl 

9v 



0+O)[lat] 
9 A 



9a 



0+[lat] 
9 a 



9°v + 



9°v 



9°a + 



Zv 



'■>A 



0.1375 4.208( 8) 
0.1390 4.492(10) 
0.1400 4.663( 9) 



3.844( 76) 
4.152(160) 
4.380(206) 



3.852( 42) 
3.978( 94) 
3.952(136) 



-1.073(49) 
-1.244(44) 
-1.150(55) 



4.925( 24) 
5.222(126) 
5.102(145) 



1.045(1) 
1.089(1) 
1.115(2) 



0.989( 19) 
1.036( 60) 
1.048(111) 



1.247(8) 
1.261(7) 
1.261(8) 



0.2530 0.2576 
0.2446 0.2491 
0.2390 0.2434 



0-(u)[lat] 

9 a 



0-(d)[lat] 
9 A 



0-[lat] 
9 A 



l-(u)[lat] 
9 A 



1 -(d) Pat] „l-[lat] 
9 A 9 A 



9° A - 



1- 
9a 



0.1375 
0.1390 
0.1400 



0.336(194) 
-0.710(251) 
0.189(257) 



-0.257(118) 
0.081(119) 
-0.129(178) 



0.592(226) 
-0.791(272) 
0.318(297) 



3.308(234) 
3.423(495) 
3.530(516) 



1.189(209) 
1.243(420) 
1.339(405) 



2.119(359) 
2.180(730) 
2.190(676) 



0.152(58) 
-0.197(68) 
0.077(72) 



0.546(093) 
0.543(182) 
0.533(165) 



previous lattice studies, the axial charge of the positive 
parity nucleon shows little quark-mass dependence, and 
they lie around the experimental value 1.26. 

We finally show the axial charges of the negative-parity 
nucleon resonances in the lower panel in Fig. [3] One finds 
at a glance that the axial charge g ^ of AT* (1535) takes 
quite small value, as g^ ~ O(0. 1) and that even the 
sign is quark-mass dependent. While the wavy behav- 
ior might come from the sensitiveness of g r~ to quark 
masses, this behavior may indicate that is rather 
consistent with zero. These small values are not the 
consequence of the cancellation between u- and d-quark 
contributions. The u- and d-quark contributions to g^~ 
are in fact individually small, which one can find in the 
columns named as g\ (")P at l anc l g° A ( d )[ lat l m Table HI 
We additionally make some trials with lighter u- and d- 
quark masses at k=0.1410. Since we have less gauge con- 
figurations and the statistical fluctuation is larger at this 
kappa, we fail to find a clear plateau in the effective mass 
plots of the two-point correlators and the extracted mass 
E± of the negative-parity state cannot be reliable. Leav- 
ing aside these failures, we try to extract g^~ . The result 
is added in the lower panel in Fig. [3] as a faint-colored 
symbol, which is consistent with those obtained at other 
k's. On the other hand, the axial charge gj[~ of N* (1650) 
is found to be about 0.55, which has almost no quark- 
mass dependence. The striking feature is that these axial 



charges, g^~ ~ and g^ ~ 0.55, are consistent with the 
naive nonrelativistic quark model calculations [HI, [Tij ]. 
g®^~ = — i and g 1 ^ = |. Such values are obtained 
if we assume that the wave functions of iV*(1535) and 
iV*(1650) are \l = l,S=\) and \l = 1,5= §) neglect- 
ing the possible state mixing. (Here, I denotes the orbital 
angular momentum and S the total spin.) 

In the chiral doublet model [j], Q , the small g^ N is 
realized when the system is decoupled from the chiral 
condensate (t/>^>). The small g^~ of TV* (1535) then does 
not contradict with the possible and attempting scenario, 
the chiral restoration scenario in excited hadrons [2|. If 
this scenario is the case, the origin of mass of TV* (1535) 
(or excited nucleons) is essentially different from that of 
the positive-parity ground-state nucleon 7V(940), which 
mainly arises from the spontaneous chiral symmetry 
breaking. However, the non-vanishing axial charge of 
iV*(1650) unfortunately gives rise to doubts about the 
scenario. 

In order to reveal the realistic chiral structure, studies 
with much lighter u,d quarks will be indispensable. A 
study of the axial charge of Roper, as well as the inclusion 
of strange sea quarks could also cast light on the low- 
energy chiral structure of baryons and the origin of mass. 

Conclusions. In conclusion, we have performed the 
first lattice QCD study of the axial charge g^ N of 
N* (1535) and A^*(1650), with two flavors of dynam- 
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ical quarks employing the renormalization-group im- 
proved gauge action at (3=1.95 and the mean-field im- 
proved clover quark action with the hopping parameters, 
k=0.1375, 0.1390 and 0.1400. We have found the small 
axial charge of AT* (1535), whose absolute value seems 
less than 0.2 and which is almost independent of quark 
mass, whereas the axial charge of Y*(1650) is found 
to be about 0.55. These values are consistent with the 
naive nonrelativistic quark model predictions, and could 
not be the favorable evidences for the chiral restoration 
scenario in (low-lying) excited hadrons. Further investi- 
gations on the axial charges of N* (1535) or other excited 
baryons will cast light on the chiral structure of the low- 
energy hadron dynamics and on where hadronic masses 
come from. 

All the numerical calculations were performed on NEC 
SX-8R at RCNP and CMC, Osaka University, on SX-8 at 
YITP, Kyoto University, and on BlueGene at KEK. The 
unquenched gauge configurations employed in our anal- 
ysis were all generated by CP-PACS collaboration 0. 
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